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Abstract 

in) 

In this paper we study the form factors for the half-BPS operators 0\ and 
the M = 4 stress tensor supermultiplet current T AB up to the second order of 
perturbation theory and for the Konishi operator K, at first order of perturbation 
theory in the M = 4 SYM theory at weak coupling. For all the objects we observe 
the exponentiation of the IR divergences with two anomalous dimensions: the cusp 
anomalous dimension and the collinear anomalous dimension. For the IR finite 
parts we obtain a similar situation as for the gluon scattering amplitudes, namely, 
apart from the case of T AB and K, the finite part has some remainder function 
which we calculate up to the second order. It involves the generalized Goncharov 
poly logarithms of several variables. All the answers are expressed in terms of the 
integrals related to the dual conformal invariant ones which might be a signal of 
integrable structure standing behind the form factors. 
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1 Introduction 

Much attention in the past few years has been dedicated to the study of the planar limit 
for the scattering amplitudes in the = 4 SYM theory. It is believed that the hidden 
symmetries responsible for integrability properties of M = 4 SYM completely fix the 
structure of the amplitudes (the S-matrix of the theory) [TJ [2]. One of the possible views 
on this subject is that the answers for the amplitudes are expressed in terms of the scalar 
integrals which are pseudo-conformal invariant in momentum space [3] which appear in 
the unitarity-based calculation of the scattering amplitudes pioneered in papers jl]. 

The (dual) conformal symmetry at weak coupling regime can be extended to the M = 
4 supersymmetric version and can be fused with the original M = 4 superconformal 
symmetry to the so-called Yangian symmetry [S] which is governed by Yangian infinite 
dimensional algebra. The Yangian like symmetries are common features of the integrable 
systems [I]. 
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At strong coupling the computation of the amplitudes in Af = 4 SYM can be reduced 
via AdS/CFT to the computation of the open string scattering amplitudes in AdS§, with 
strings ending on D3-brane positioned at some fixed value of the radial AdS§ coordinate z, 
in the quasi-classical regime [6] which in turn can be formulated as the problem of finding 
the minimal surface in AdS§ with special boundary condition (see [7] for review). This 
problem has recently been reduced to that of solving the set of functional equations for 
the conformal invariant cross ratios as functions of the spectral parameters- the so-called 
Y-system [5|. The Y-systems usually appear in integrable systems [H] which is another 
hint that the Af = 4 amplitudes have some underling integrable structure. 

In strong coupling regime the natural generalization of the Y-system for the amplitudes 
is the Y-system for the form factors [10]: the matrix elements of the form 

(0|O|^...^>. (1.1) 

where O is some gauge invariant operator which acts on vacuum and produces some state 
l-Pi 1 • • -Pn n ) w ith momenta pi . . .p n and helicities Ai . . . A„. In the dual string theory this 
matrix element can be described via the open string scattering amplitudes, with strings 
ending on D3-brane positioned at some fixed value of the radial AdS§ coordinate z, in 
the presence of some closed string state [10] • 

One can wonder whether these objects at weak coupling possess similar features as 
the amplitudes or in other words whether form factors are influenced by the Yangian 
symmetry (or some analog of it) and whether they are fixed by it. Also the general 
structure of the form factors at weak coupling should be understood. 

Being inspired by the two-loop calculation of the form factor associated with the 
operator Vx from the stress-tensor super conformal multiplet of the Af = 4 SYM theory 
performed long ago by van Neerven [IT] we would like to study systematically some types 
of form factors in planar Af = 4 SYM at weak coupling for half-BPS operators Oj and 
the Konishi operator /C. For the former type of operators there was recently an interest in 
studying the correlation functions and their connection to the amplitudes and the Wilson 
loops [12]. A new kind of relation has been proposed between the logarithm for such 
correlation functions and the double logarithm of the MHV gluon scattering amplitudes. 

The dual-conformal symmetry plays an important role in another remarkable property 
of the Af = 4 SYM - the Wilson loop/ Amplitudes duality [T3l_T4]. In this duality (we re- 
strict ourselves to the most studied case of the MHV amplitude! sector) the dual-conformal 
symmetry is understood as conformal symmetry of light-like Wilson loop constructed of 
the segments which satisfy the following property: 

X i,i+1 = X i ~ X i+1 = Pi 1 (1-2) 



1 MHV (maximally helicity violating) amplitudes by definition are called the amplitudes with all par- 
ticles being treated as outgoing and the net helicity Ag being equal to n — 4 where n is the number of 
particles. 
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where pf are external momenta of the dual MHV amplitude. The dual-conformal symme- 
try is broken on-shell for the amplitudes due to the presence of the IR divergences (these 
IR divergences correspond to the UV divergences for dual Wilson loops); however, the vi- 
olation of the dual-conformal symmetry is controlled by the 1-loop exact anomaly, which 
in turn can be used to make constraints for the finite part of the corresponding MHV 
amplitude, i.e. one can write the anomalous Ward identities allowing one to constraint 
the finite parts (see, for example, review [7J for details): 

n n g.2 

Y,(2<Xidi - x^)Fm[log W n ] = T cusp £ log -J**-^, (1.3) 

8=1 i=l X i-l,i+l 

where T cusp is the so-called cusp anomalous dimension [T5] known from a solution of 
the integral equation [16]. These identities can fix the finite part for a small number of 
legs/cusps (namely, one can fix it for Wilson loops with n = 4 and n = 5 cusps) [7J. 
In fact, the famous BDS conjecture was the simplest possible ansatz of these identities, 
which is not precisely correct for a number of external legs > 5. 

One may wonder if there is a similar duality for the form factors/ Wilson loops and one 
can use similar arguments to obtain information on the finite parts of the form factors. 
We hope that our calculation shades some light on the possibility of such duality. 

The paper is organized as follows. In Sect. El we present the general considerations of 
the form factors in the M = 4 SYM theory and introduce the operators to be discussed 
later. In Sect. EJ we study the form factors for both protected and non-protected operators 
with naive conformal dimension A = 2 and confirm the results obtained long time ago in 
In Sect. IU we study the half-BPS operators Oj for arbitrary conformal dimension 
Aq = n and present the one- and two-loop calculations of the corresponding form factors 
which suggest the exponentiation of the IR divergences. Also, in the same section, we 
discuss the collinear limit for which the finite parts take a simple form. In Sect. we 
discuss in more detail the dual conformal invariance of the integrals contributing to the 
calculation of the form factors. We conclude with some remarks concerning the form 
factors and Wilson loop duality. In the appendices we give the details of our calculations. 
Appendix. A contains the Lagrangian of the M = 4 SYM theory together with the 
Feynman rules. In App. B we present the analytic expressions for the integrals entering 
into our calculations with their e-expansion. The results of our work have been reported 
at the international conference devoted to the memory of A.N. Vasiliev held in Saint- 
Petersburg on 18-21 October 2010. 

2 General considerations 
2.1 Form factors in M = 4 SYM 

Consider the Lagrangian £jv=4(W) for the Af = 4 SYM theory coupled to some external 
classical current J through some gauge invariant local operator 0[VV] (for all the details 
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concerning the explicit expression for the Lagrangian together with Feynman rules we 
refer to App. A) 

C M=4 (W) C M=4 (W) + 0[W]J, (2.4) 

where we collectively refer to the whole M = 4 on-shell multiplet as W which consists of 
the physical gluon states with positive and negative helicities, four gauginos A^ with 
positive and negative helicities and also six real scalar states <Pnmi where N and M are 
the SU(4)r indices, which can also be re-arranged into 3 complex cj) 1 scalars, I is an index 
of SU(3) subgroup of SU(4:) R . By default we assume everywhere the planar limit. 

Then one can study the following processes where the operator O acts on the vacuum 
and produces some state \p^ . . -Pn n ) with momenta p\ . . .p n and helicities Ai . . . A n 

(0\O\p^...p x n »). (2.5) 

Schematically, it is shown in Figfj] 

q 



Pi P 2 Pn-1 Pn 

Figure 1: Feynman diagram for the matrix element of the operator O 

This is a general situation in QFT and one can keep in mind, for example, 7* — > Jet's 
process [TTJ where we take into account all orders in a s but the first order in a em . In 
perturbation theory the latter type of processes can be thought of as the matrix elements 
of the following form: 

(0b e Q m CD |^...^), (2.6) 

where j2^ D is the QCD quark electromagnetic current. 

The matrix element (OlOlpi 1 . . -Pn n ) m some sense can also be viewed as the general- 
ization of the scattering amplitudes, which in " all ingoing" notation can schematically be 
written as (O^ 1 . . -Pn n )- 

In the language of the dual string theory, in the M = 4 SYM case this process can 
be described as an insertion of some close string state (which corresponds to O local 
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operator) on the worldsheet in addition to n open string states (which corresponds to 
IPi 1 ■ ■ -Pn n ) state in the dual theory). 

For the construction of particular examples of the objects discussed above we choose 
the following set of the gauge invariant operators (we use the component notation of the 
M = 4 SYM), the lowest stress tensor supermultiplet members: 

Cjj = Tr(0 7 0j), 

Vj = Tr(0 J 7 ), (2.7) 

with naive mass dimension Ao = 2 which coincides with the conformal dimension due to 
the lack of quantum corrections. These operators can be viewed as the lowest members 
of the stress-tensor multiplet 

T AB = Tr (w A W B - ^5 AB W c Wc^j , (2.8) 

where A, B, . . . = 1, . . . , 6 are the SO(6)r ~ SU(4:) R indices, I, J, . . . = 1, 2, 3 are the 
indices of SU(3) subgroup of S't/(4)^, and W A is some constrained chiral superfield in 
M = 4 superspace containing all components of the M = 4 supermultiplet. 
Other objects are the so-called half-BPS operators 

Of = Tr(^), (2.9) 

whose naive mass dimension coincides with conformal dimension A = n being protected 
from the quantum corrections, and the lowest component of the Konishi supermultiplet 

/C=]TTr(0V/), (2.10) 

7=1 

with naive mass dimension Ao = 2 and has nonvanishing anomalous dimension due to the 
presence of the UV divergences. The calculation of this anomalous dimension has been 
intensively discussed during the last few years [HI [19] . 

Since the Konishi operator is not protected, the corresponding form factors a priori 
do contain the UV divergences and hence must be UV renormalized. It means that one 
has to consider the renormalized form factor 

(0\}C R \p^...pt), (2.11) 

where 

K-R - z k K-b- (2.12) 

Here is the renormalization constant which appears due to the UV divergences and 
which should be calculated to the same order of perturbation theory as the form factors. 
After such UV renormalization we are left only with the IR divergences. All the statements 
concerning the Konishi operator are valid for the renormalized one. 
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We choose for simplicity the state \p^ . . -Pn") produced by the operator O to consist 
of scalars only and then we can write it as \pi ■ ■ ■ p n ) without helicities. We also restrict 
ourselves to the states with the number of particles equal to the naive mass dimension of 
the operator O, i.e. we consider the states consisting of Ao scalars. 

2.2 Calculation strategy 

For the calculation it is convenient to use the M = 1 formulation of M = 4 SYM and 
perform an explicit computation in terms of the Af = 1 superfields in momentum space. 
Our computation is familiar, from a diagrammatic point of view, to perturbative compu- 
tations of anomalous dimensions [15] . However, there is a significant difference: each of 
our supergraphs is UV finite except for one. So all the divergences that appear throughout 
the calculation have the IR nature. 

The operators O = {C u , V/, /C, C»} n) } can be considered as the lowest components of 
the following M = 1 local operators: 

C u = Tr($i$j),I?J 

V/ = Tr(e- 9 W$/),/^J 

Oj ] = Tr($?), 

K = £T!r(e- ffV $V v $j), (2.13) 

where are chiral N — 1 superfields, and V is M = 1 real vector superfield (see App. A 
for details). The operators Cjj, Oj^ are chiral and V/, K, are non-chiral from the N = 1 
supersymmetric point of view. 

We use the following notation for the form factor of the corresponding operator 

F(p 1 ...p n ) = {p 1 ...p n \O(q)\0). (2.14) 

We expect the following factorization property for T to hold: 

?{P\ ...p n )= FtreeiPl ■ ■ -Pn){^ + loops), (2.15) 

where J-'tree stands for the tree level contribution, and "loops" schematically denote the 
contributions of the next orders of PT. It is convenient to consider the ratio 

M = = (1 + loops) = \ l M {l \ 

tree \—q 

where A = g 2 N c is the 't Hooft coupling which stays fixed when N c — > oo. 

Consider first the chiral case. To calculate the form factor it is convenient to consider 
the generating functional for the one-particle irreducible super diagrams r[$ c/ , J] in M = 1 
superspace. It can be obtained from the generating functional 

Z[j,J] = J V^ I ,V,...)exp[S Af=i + J d e zJ(z)0{z) + J d G zTr{j{z)$(z))), 
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after Legendre transformation with respect to external chiral sources j (note that the 
source J is untouched). After performing the D-algebra each supergraph gives a local 
contribution in #'s, and r[$ d , J] can be written as ( we imply the mass shell condition 
pf = when performing the D-algebra) 



J] =^A / r w [$ d , J] (2.16) 
1=0 

= E A ' / d 'Pi ■ ■ ■ d 'Pn d 6 z J(-q, 9)Tr 0) . . . $ d (-p„, 0)) .M (Z) (pi, . . .p n ) + 0(J 2 ), 

1=0 

where d 6 z = d 4 qd 2 9, M."' is given by the sum of scalar integrals. Thus, 

^n+lp(0 

We stress that on-shell condition pf = and momenta conservation q + pi + . . . + p n = 
are implemented to obtain the latter expression. 

The situation is a bit more involved in the nonchiral case. All the integrals in 
r[$ d , $ d , J\ {J is a non-chiral source) are now in full M = 1 superspace / d 8 z, where 
d 8 z = d A qd 4: 8 and the expression for T[$ d , $ d , J\ contains extra terms 

r[$ d , $ d , j] = £A'r ( °[$ c ',<j> d , j] = (2.i8) 

1=0 

= E A' / • • • ^ 4 Pn ^(-?, 0, 0) [Tr (<l d (-p 1 , 0) . . . $ d (-p„, 0)) A^^Cpx, ...Pn) 

1=0 J 

+Tr (^^(-px, 9)... D a &\-p n , 9)) Mf a (p u ...p n ) 

+Tr (D 2 $ d (-p u 9)... D 2 $ d (-p n , 9)) Mf{p u . . . Pn)} + 0(J 2 ). 

From the point of view of Af = 1 superspace the additional terms correspond to 
the operators of higher dimension and one actually has a mixing of several operators. 
However, from the point of view of components, one can always consider a projection on 
a particular component of a superfield and we choose the scalar component insofar. Then, 
the last terms of eq. (12.181) are irrelevant for our calculation and can be dropped. 

We perform all the calculations in the formalism of M = 1 superspace and at the end 
take the projection to 9 = 9 = 0. There are pluses and minuses of this approach. The big 
advantage is the drastic reduction of the number of diagrams compared to the component 
case together with the simplified form of the scalar integrals. Its disadvantage is that 
we do not use the power of the on-shell M = 4 covariant methods used in perturbative 
studies of the amplitudes [20] (see also recent [21]). The application of this method for the 
calculation of the form factors when some legs are off-shell requires some modification. 
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2.3 IR finite observables based on form factors 

As the amplitudes, the form factors (0| Olp^ 1 . . -p^ n ) with on-shell momenta are, strictly 
speaking, ill-defined in D = 4-dimensional space-time due to the presence of the IR 
divergences, and, hence, some IR regulator must be introduced - in our case it is the 
parameter fi coming from the dimensional regularization which also breaks the conformal 
symmetry. In other words, one may say that (0| C^l^^ 1 . . . p^ n ) are the intermediate objects, 
and the true physical quantities are the IR safe observables constructed of (OlClp^ 1 . . . p* n ) 
and which are free from the IR regulator (see, for example, the discussion of the IR finite 
observables for M = 4 SYM and M = 8 SUGRA in [221 [23]). Indeed, as in QCD for 
7* — > Jet's processes we are really interested in the total cross section a to t{l* — > Jet's) 
or some differential distributions rather than in the matrix elements (Olj^^lpi 1 . . -p^) 
themselves. This kind of observables are IR finite due to the Kinoshita-Lee-Nauenberg 
(KLN) theorem which states that it is not sufficient to consider only the processes with 
the fixed number of final particles. To get the physical result, one has to include all 
the processes allowed by conservation laws in the same order of perturbation theory 
with emission of extra soft quanta and integrate over their momenta. Practically, if the 
dimensional regularization is used (the IR divergences manifest themselves through the 
appearance of the 1/e poles), part of the poles cancel between the loop integrals from the 
virtual contributions and the phase space integrals from the real contributions coming 
from the processes with additional particles, while others are absorbed in the functions 
describing probability distributions of the initial and final states (in [221 123] we call them 
initial-collinear and final-collinear divergences which appear as a collinear configuration 
of initial and final particles). 

Consider, for instance, the total cross section a to t for the process 

J — > anything from Af = 4 supermultiplet 

for classical current J coupled to M = 4 through some local gauge invariant operator O. 
Due to the optical theorem 



o-tot(s) ~ -Im s 
s 



d n xexp(-iqx)(O(x)O(0)} 



, q 2 = -s, (2.19) 



The two-point function for the operators O apart from the canonical mass dimension A 
can have anomalous dimension 7 = 7(A) being a function of the coupling constant 

^)Om ~ (x2)(Ao 1 (1 _ £)+7r (2.20) 
After some calculation this gives the total cross-section 

atot{s) ~ r(A + 7)r(A + 7 -i) s 3-Io-,' ( 2 - 21 ) 
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and its asymptotic at weak and strong couplings can be studied (compare this with C.3 
from [21]). 

In M = 4 SYM, as in any conformal theory, the operator O is protected, which 
means that it does not receive quantum corrections and 7 = 0. Then the cross section 
is independent of the coupling constant and behaves like ~ C/s 3 ~ A °. From the latter 
expression it might seem that we get violation of unitarity since we can get increasing 
cross sections for protected operators with conformal dimension greater than 3. But it is 
not the case since the statement about the unitarity holds only for the operators which 
give rise to renormalizable interactions, i.e. with conformal dimension less than 3. 

If one is interested not in a tot but in some differential distributions, then the optical 
theorem is not very useful any more, and direct computations must be done. The form 
factors discussed here can be viewed as the building blocks in the same sense as the 
amplitudes for the inclusive cross sections. 

3 Form factors with Aq = 2 

In the following two sections we give explicit results for the direct diagrammatic compu- 
tation of the form factors of the operators introduced above in the planar limit H More 
concretely we present the results for non-chiral operators V/, JC in the leading order in 

(n) 

A and for the chiral operators Cu, 0\ in the next-to-leading order. The dimensional 
regularization (dimensional reduction to be precise) with D = 4 — 2e is used. All the 
divergences except for the specially mentioned cases have the IR (both the infrared itself 
and collinear) nature. The Feynman rules for the supergraphs are given in App. A. The 
complete list of all the necessary scalar integrals is given in App. B. 

3.1 Cjj, V/ and JC form factors at 1-loop 

The corresponding tree level and one-loop Feynman diagrams are shown in FigfSJ 

For the chiral operator Cjj, after performing the D-algebra for the supergraph CI, the 
resulting expression is 



where the scalar integral Gx(si2) 1S given in App. B Hereafter we will suppress index cl 
in $ and <J. 

2 g —> and N c —> oo so that A = g 2 N c =fixed. 




(3.22) 



where s^- = (pj +pj) 2 , so that 




(3.23) 
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CO 



C1 



B1 



B2 



Figure 2: The relevant supergraphs. The internal black lines correspond to chiral propa- 
gators wavy lines correspond to vector (V a V b ) propagator (see App. A). CO is 
the tree level diagram, and the rest are one- loop ones. External lines are $ or i>, and the 
lower bold line represents the insertion of the corresponding operator in modern notation. 
For the chiral operator Cjj only the diagrams CO and CI contribute, while for non-chiral 
operators V/ and /C the other two (Bl and B2) are also relevant. 

For the non-chiral operator V/ after performing the D-algebra one has 
CI = 2 ((-G (p 2 ) - G (pl) + G (s 12 ) + (s 12 - V \- pl)Gi(s 12 )) Tr ($ J <Z>j) 

-Gf (s 12 ) Tr {p^ J D a ^ I ) + Gi(sia) Tr (D 2 ® 3 D 2 <&^j , 
Bl = 2G (ft 2 ) Tr ($ J $/) , (3.24) 
B2 = 2 (-G (si 2 )Tr (i» J $/)+Gf (s 12 )Tr ( J D^$ J J D a $ J ) +Gi(si 2 )Tr (dWD 2 ^ , 

where all the scalar integrals are also given in App.B. So keeping only the terms that are 
relevant for our discussion, which are proportional to Tr (<I J $/) , one gets 

= Cl(s 12 ,pl,pl) + Bl(pj) + Bl(pl) + B2(s 12 ) = 2(s 12 - V \- V ^G x (s X2 \ (3.25) 

The integral G\ is UV finite which reflects the fact that Vj is a protected operator. 
For the non-chiral Konishi operator /C, after performing the D-algebra one has 



CI = 



Bl 
B2 



6 [(-Go(pl) ~ G (p 2 2 ) + G (s 12 ) + (s 12 - V \- V \)G x (s X2 )) ]T Tr 

3 3 \ 

-Gf(s 12 ) E Tr (D^¥D a ^j) + d( Sl2 ) £ Tr (d^'D 2 ^) , 

7 I J 

6G (p?) (3.26) 

7 

(3 3 

-C ( Sl2 ) £ Tr + Gf (s 12 ) E Tr (Wl)^ 
7 7 
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3 

+Gi(s 12 )^Tr (D^D 2 ®! 
i 

and again selecting the proper structures, prior to the application of the on-shell condi- 
tions, gives 

M {1) = Cl(sn,plpl) + Bl(pl) + Bl(pl) + 2B2(s 12 ) = Q(s 12 -p 2 1 -pl)G 1 (s 12 )-QG (s 12 ), 

(3.27) 

The UV divergent part of the answer is given by 6Go and extracting the coefficient of 
the 1/e pole, which is the first coefficient in the anomalous dimension expansion 7ac(A) = 
7^ A + • • • , we obtain the well-known result 

t£> = ± (3.28) 

We see that up to one loop all the form factors for the operators Cu, V/, JC, Oj are 
proportional to G±, the scalar triangle function (see App. B). 



3.2 Cu form factors at 2-loops 

We see that the form factors associated with Cu and V/ are equal to each other at the one- 
loop level. This is because Cu and Vj are different components of the J\f — 4 conserved 
stress tensor. In what follows we compute the A 2 contribution to M. for Cu since for 
the chiral operator the D-algebra is essentially simpler. The corresponding diagrams are 
shown in Fig 131 

Their contribution to the form factor are summarized in Table HJ All the relevant 
integrals are given in App.B 



N 


Cu 




CI 


2si 2 Gi(s 12 ) 


Sii+iGi(su + i) 


C2 


As\ 2 G 2 {s l2 ) 


s 2 i+l G 2 (su + i) 


C3 


2s 12 G 3 (s 12 ) + 2s l2 G A (s l2 ) 


SU+1G3 (sji+i) + Sii + iGi(su + i) 


C4 


-6s 12 G 3 (s 12 ) 


—^ s ii+lG-i{Sii + i) 


C5 







C6 


2G b 5 ( Sl2 ) 





C7 





+ Si + u +2 + Sh +2 )Gq(su + i, Sj_|_ij_|_2, Su +2 ) 


C8 





Si+li+ 2 GT(su + i, Si + u +2 , Su +2 ) 



Table 1: The contributions to the form factors from the individual diagrams 



Adding all together and combining with the leading order one gets 
M {2) = C2 + C3 + CA + C5 + C6 = 2G% + 2G\ - 4si 2 G 3 + 2s X2 G A + As\ 2 G 2 . (3.29) 

12 




C5 C6 C7 



Figure 3: The relevant supergraphs in the chiral case. CI is the one-loop diagram, and 
the rest are two-loop ones. For the chiral operator Cu with two legs the last two diagrams 
CI and C8 do not exist, they are only relevant for the operator O n with n > 3. A grey 
circle is the one-loop effective vertex. 

Using the identity 

b s u 
^5 + ^5 = 2S12C3 — S12G4 + 

one can reduce it to As\ 2 G2 + s\ 2 G5 and finally get 

M = 1 + A(2s 12 G x ) + \ 2 (As\ 2 G 2 + s\ 2 G 5 ), +0(A 3 ). (3.30) 

The general structure of form factors Ai "with two external legs" in the gauge theory 
with zero beta-function has the following forrrH: 

log(-M) = \il (M(s t , t+1 /fi 2 )) + 0(e). (3.31) 
A i=i 

Here we introduced 




(3.32) 



where r£2 are the coefficients of perturbative expansion of the cusp anomalous dimen- 
sion T cusp (X) = J2i r^spA' which is a universal quantity that governs the IR behavior of 

3 The IR exponentiation for two-leg form factors in QCD was established earlier in |25j . 
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gauge theory amplitudes and the UV behavior of the Wilson loops, and some local gauge 
invariant operators, are the coefficients of perturbative expansion of the so-called 
collinear anomalous dimension G(X) = Y^iG^X 1 and are some constants. The quan- 
tities G"> and C"> are regularization and scheme dependent. Performing the expansion 
of the integrals G±, G2, G5 in e (see App. B) and introducing the notation 

a = ^e-, (3.33) 
where je is the Euler-Mascheroni constant, one gets 

log(M) = a (£H) £ + C2 ) + J M (| + |) + 0(0 3, ,3,4) 



where ( n are the Riemannian zeta functions 

00 



00 

fe=l 



From this answer and comparing with eq. fl3.32p we can extract the first two terms of 
perturbative expansion over a for the cusp and the collinear anomalous dimensions and 
the finite terms 

= 4, rf u ) p = -8Ca, (3.35) 
= 0, G (2) = -C 3 , (3.36) 
C« = -C 2 , C (2) = 0. (3.37) 

Note that the maximal transcendentality principle holds which in our case means 
that if we attach to each logarithm and 7r the level of transcendentality equal to 1 and 
to polylogarithms Li n (x) and ( n the level of transcendentality equal to n, then at the 
given order of perturbation theory the coefficient for the n-th pole l/e n has the overall 
transcendentality equal to 21 — n, where I is the number of loops. For a product of several 
factors it is given by the sum of transcendentalities of each factor. 

The leading IR behavior of M. in this case can also be captured by considering the 
Wilson line with one cusp [26] . So in this sense the dual description in terms of Wilson 
loops for such form factors is well known. 

One can see that the finite part for the form factor is given only in one loop and 
vanishes at two loops. However, this is a scheme dependent result, and, for example, if we 
choose a different scheme and replace exp(/e7^) for the Z-th loop by T(l — e) 1 , we obtain 
in this scheme: 

(7« = 0, C {2) = -Cl (3.38) 

while the first two coefficients in the perturbation theory for the cusp anomalous dimension 
^iusp an d ^cusp remain the same, which reflects the fact that they are scheme independent. 

The same result is true [TT] for the finite part for the form factor of a slightly different 
operator Vx but belonging to the same stress-tensor superconformal multiplet. 
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4 Form factors with Aq = n, n > 2 



CO 

Figure 4: The tree contribution to Oj. 

Here we present the results of calculation of the form factors of the chiral half-BPS 
operators Oj introduced earlier. The tree-level contribution for the form factor is pre- 
sented on Fig. HJ In the first order of perturbation theory, similar to the form factors 
of operators with conformal dimension 2, the contribution is given by the triangle type 
diagram and the corresponding form factor, after performing the D-algebra and the color 
algebra, is 

n 

■M (1) = £s M+ iGi, ( 4 - 39 ) 

1=1 

where we assume hereafter s n+ i <n+ i + i = s^+i. 
(n) 

4.1 Oj , n = 3 form factors at 2-loops 

At the second order of perturbation theory the corresponding diagrams are shown in Figj3] 
and their contributions are summarized in Table [TJ 

n 

J\y[( 2 ) — (s 2 ii+1 G2(s ii+ i) — Sjj +1 G 3 (Sjj +1 ) + Sn + iG4(Sii + i)^j 

i=l 

n 

+ ( S i+U+2G7(Sn + i, Si + u + 2, Sa + 2) + Sjj +1 G 7 (,Sj +1 j + 2, Sa + i, Sn + 2)) 

i=l 
n 

Sii+1 + Si+li+2 + Sii+2) Gq(s 

) ( 4 - 4 0) 

i=i 

The next step is to establish the factorization property (I3.31|l3.32p . Expanding the 
relevant scalar integrals in e (see App. B), we obtain for log(A'l): 

log(M) = ± a ^y^ + 1) + g a2 (to) (6 + Z|) +fln . part . (4 . 41 ) 
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As in the case of the form factors of the operators with conformal dimension 2, we can 
extract the first two terms for the cusp and collinear anomalous dimensions. This gives 

r« p = 4, rfj sp = -8C2, (4.42) 

G« = 0,G (2) = -7C 3 . (4.43) 

Notice that the values of the cusp anomalous dimension are universal and coincide 
with f)3.35p . while those of the collinear anomalous dimension depend on the form factor 
at hand. 

We would like to emphasize the highly nontrivial cancelations between the polylog- 
arithms that occurred for log(A^) for the whole set of scalar integrals (the individual 
contributions to the poles from the scalar integrals are usually complicated polynomials 
of logarithms and polylogarithms of different weight, see, for example, the expansions in 
e of G G and G 7 in App. B). 

We see that the IR factorization property holds for the form factors like for the am- 
plitudes. 

(n) 

4.2 Oj form factors for n > 3 at 2-loops 

The corresponding contribution to the form factor up to A 2 is similar to the case of n = 3 
but has an additional term coming from the factorized diagrams 

n 

j^m _ (s 2 i+1 G 2 (si i+ i) - s u+ iG 3 (s ii+1 ) + s <i+1 G 4 (sii+i)) 
1=1 

n 

+ ( S i+li+2G7(Sa + i, Sj + ij_|_2, Sjj+2) + Su + iGj(Si + ii + 2, Sfi+l, Sii+2)) 

i=l 
n 

Sii+l + Si+li+2 + Sii+2) Gq(Sh + \, Sj+ij+2, Sii+2) 

i=l 

n n 

+ ^2^2su +1 G 1 (s u+ . 1 )s jj+ iG 1 (s jj+ i) (4.44) 
i=i j=i 

Performing the expansion over e we obtain the logarithm of the form factor up to the 
second order of perturbation theory \og(Ai) 



log{M) = g a ^ (_ 1 (^)" 2< + f j + Fi„.pa rt . 

(4.45) 

The first two coefficients for the cusp and collinear anomalous dimension which we 
can extract from the above expression coincide with the coefficients obtained earlier for 
n = 3, (I4.42p and f)4.43p . respectively. As for the finite part 

Fin.part. = AF«(s 12 , . . . , s nl ) + \ 2 F^(s 12 , • • • , s nl ) + 0(A 3 ), (4.46) 



16 



at one loop it is trivial F^ 1 ' = 0, and the two loop expression F^ 2 \ contrary to the previous 
case, is a complicated function containing logarithms, polylogarithms and generalized 
Goncharov polylogarithms [30J of several variables. All the relevant expressions can be 
found in App. B. 

Note, the result is still much simpler than in the non-supersymmetric case [29] and 
the maximal transcendentality principle still holds. 

4.3 Collinear Limit 

Here we restrict ourselves to the three-leg form factors for which we can study the simpli- 
fied kinematics and express the finite part in terms of logarithms only without polyloga- 
rithms or Goncharov generalized polylogarithms. The most difficult part which appears in 
our calculation comes from the diagram involving the interaction of three external fields. 
It reduces to the integral Gj which is expressed in terms of the Appell function of two 
variables 

F 1 (l;2e,l;2 + e\x,y) 

and after the e-expansion one obtains the generalized Goncharov polylogarithms [3D]- One 
can see from integral representation of the Appell function 

H*'MMx,y) = nanc _ a) J o {1 _ ux)bi{1 _ uy)b2 du,Rea,Re(e-a)>0 (4.47) 

that the only way to achieve the desired simplification is to have one of the arguments 
equal to or to 1. In the first case one gets 

Fi(a;6i,6 2 ;c|x,0) = 2 Fi(o, b x ; c\x), (4.48) 

and similarly in the second case 

F 1 (a;b 1 ,b 2 ;c\x,l) = r(c)r(c - a - b 2 ) g _ 

T(c - a)T(c - b 2 ) 

Such a simplification can occur in two-dimensional kinematics when one of the kine- 
matical variables Si2,si3 or S23 equals 0. The other motivation for this kinematics is the 
recent strong coupling calculations which have been performed for the AdS$ sub-manifold 
of AdSz, which corresponds to the degenerate 1 + 1 kinematics in a dual theory [10J. 

The 1 + 1 dimensional kinematics necessarily contains a collinear configuration of the 
space components Pi of momenta pj. For massless gauge theory it is known that in such 
collinear limit the factorization of the IR divergencies fails. For the partial color ordered 
amplitudes in collinear limit when two momentums pi and pi + \ are replaced by zp and 
(1 — z)p the deviation from the factorized form is governed by the so-called "loop splitting 
functions" r)P(e, z,p 2 ), I being the number of loops. In the J\f = 4 SYM theory r®(e, z,p 2 ) 
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have an iterative structure, so one can write the following relation valid in collinear limit 
(see, for example, the discussion in [15] ) 

log(M„) -> ijWU + £ \ l T%rf(lt, z,p 2 ) + £ \'F«\ C °" + 0(e) 



e 2 \/i / \ sin(7re) \ 2 / ^ 1 — z 



We expect that similar violation of the IR factorization happens in the case of the 
form factors. Indeed, in the S23 — > limit we have, up to A 2 



MM) = £ a - (-1 + |) + £ a> (te ) (A + I j (4 , 0) 



5 Dual conformal invar iance 

Here we would like to discuss the property of dual conformal invariance of the integrals 
appearing in our calculation. It is believed that all the integrals entering into the calcula- 
tion of the amplitudes should be dual conformal invariant when external legs are off-shell 
[3] . This means that the answers should be expressed in terms of the conformally invari- 
ant cross-ratios, which restricts the form of the result. On mass-shell this dual conformal 
symmetry has an anomaly but still remains a very important ingredient in understanding 
the properties of the on-shell amplitudes and the Wilson loops (see, for example [21 E])- 

When calculating the form factors one has similar integrals though they do not possess 
explicit dual conformal invariance. However, it is remarkable that all the integrals that 
appear in our calculation can be obtained from those contributing to the amplitudes by 
some reduction which we describe below. Because of this connection they also bear some 
conformal properties. 

Consider several examples. At one loop there is a single triangle diagram contributing 
to all the form factors. The one-loop triangle is the first in a chain of the ladder type 
diagrams [31] and has the property of dual conformal invariance [311 [32]. This diagram 
is to be connected to the box diagram, which is dual conformal [32]. in the following way. 
Consider the one-loop off-shell box diagram in momentum space which is given by the 
integral 

f d D k 1 

D 1 - loop (p 1 ,p 2 ,p 3 ,p 4 ) = / s Du2(u ^77— wTT — -V2- (5-51) 

J {2ti) u k 2 (k - pi) 2 {k + p 2 ) (k + P2+ P3) 
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Introducing the dual coordinates 

Pi = %12, P2 = £23, P3 = x 34, Pa = X41, k = x 5 
we can rewrite the initial integral in the following form 

d D x 5 1 



D 1 - Lo °p(x 1 ,x 2 ,x 3 ,x,) = J 



2222 

<-y» iJ nn" /yi" /y>" 

x 15 x 25 x 35 x 45 



,-y»2 ,-y»2 

X 13 X 2 4 



(5.52) 



where we introduced the notation and F) is the function given in 

X and Y are the conformal cross-ratios 



X 



x 12 x 34 

2 2 ' 
x 13 x 24 



F 



,7^2 ™2 
x 14 x 23 

/-y»2 . 2 

X 13 J-24 



If we now multiply (I5.52p by x\ 2 and take the limit x 2 — > 00, we obtain the one-loop 
triangle diagram [32] 



C 1-/o OP = Hm ^2 



X2— ><X) 



12 



2 2 2 2 

iy nn" rf" nn" 

x 15 x 25 x 35 x 45 



<i 4 £5 



2 2 2 
x 15 x 35 x 45 



x 34 



(5.53) 



with 



X 



x -f,y 

x 13 



/-y>2 

x 14 
x 13 



Schematically, the described procedure of obtaining the one-loop triangle diagram 
from the one-loop box diagram is represented in Fig. [5j On the left hand side one 



1 

A 

/ V 

/ v 
/ \ 
/ \ 

5 \ 



Figure 5: The one-loop triangle diagram from the one-loop box diagram. The red dot 
should be taken to infinity, and the blue line (propagator in momentum space) should be 
contracted to a point. 

has the one-loop box diagram together with the dual grid, the black lines represent the 
denominator terms appearing in the integral in x-space. Taking the limit x 2 — > 00 in 
(15.53P is equivalent to removing the grid line £25 from the dual graph and shrinking the 
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crossed line to a point in the initial graph. The resulting initial graph corresponds to the 
triangle diagram, as is shown on the right hand side. This way the triangle diagram can 
be obtained from the box one and inherit its property of dual conformal invariance. 

In the same manner one can show how the other diagrams which appear in our cal- 
culation can be obtained from the corresponding diagrams entering into the amplitude 
calculations. Schematically, we present this procedure in Figj6j 





Figure 6: The two- loop ladder type triangle diagrams from the two- loop box diagrams. 
Green arc corresponds to the presence of a numerator 



For the vertical box diagram one has to take the limit £3 — > 00. As in the previous case 
this corresponds to removing the grid line X36 (and shrinking the corresponding crossed 
line) which results in the diagram shown on the right hand side. This is exactly the ladder 
integral that appears in two-loop calculation of the form factor with n > 2 legs. 
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For the horizontal box diagram one should take the combined limit X2, x% — > oo which 
is schematically shown on the right hand side. This is the new type of integrals which 
appears only in the case when n > 2. 

The same procedure is expected to work at higher levels of perturbation theory. Our 
conjecture is that the integrals appearing at any order of perturbation theory are obtained 
from dual conformal invariant diagrams by contraction of n propagators at the n-th loop 
order. 

6 Discussion 

In this paper we continue the perturbative study of the form factors at weak coupling 
for the M = 4 SYM theory which was initiated in [TT] where the author considered the 
form factor for the operator Vx of conformal dimension 2 which in M = 1 superspace 
notations takes the form 2Tr$i$i — Tr$2$2 — Tr$3$3. The original calculation has been 
performed in components and the form factor was computed up to the second order of 
perturbation theory. In our paper, we started with the operators belonging to the stress- 
tensor super conformal multiplet, namely, with V/ = Tr y<j) J ' (p^ and Cjj = Tr (0/0j) and 
calculated them up to the first and second order of perturbation theory, respectively. We 
obtained the same results as in [11] as it was expected. 

Then we considered the Konishi operator /C = J2i Tr (0 7 0jj with classical conformal 
dimension 2 in the one loop approximation. Not being protected by supersymmetry this 
operator has the UV divergences which have to be renormalized. 

The main result of our paper is the calculation of the two-loop form factors for the 
half-BPS operators 0\ n ^ = Tr (0™) , n > 2. At the one loop level the answer for the form 
factor is very simple given by triangle diagram while at two-loops it is essentially more 
complicated. The analytical expressions for the two-loop results are given in terms of the 
Gauss hypergeometric functions and the Appell function of two variables. Their expansion 
over e up to 0(e) leads to logarithms, polylogarithms and, because of the Appell function, 
generalized Goncharov polylogarithms of several variables. However, all of them have the 
same transcendentality [16j [33] . 

In the simplified kinematics the answers become much more simple. Thus, in two- 
dimensional (or 1 + 1-dimensional) kinematics for the form factors of the half-BPS oper- 
ators 0\ it is possible to get rid of the Appell functions and after expanding over e to 
get the result in terms of the ordinary logarithms. 

For all the considered form factors we observe the factorization of the IR divergences up 
to the second order of perturbation theory. This allows us to derive the first two terms of 
expansion for the cusp anomalous dimension in coincidence with the other calculations and 
for the collinear anomalous dimension, where we obtained the first nontrivial coefficient 
at two loops = —7(3. It differs from the collinear anomalous dimension coming 

from the amplitude calculation but coincides with collinear anomalous dimension for the 
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light-like Wilson loop [3 |26]. 

The remarkable part of our calculation besides factorization is the fact that the one- 
and two-loop integrals contributing to the form factors of the operators 0\ n \n > 2 
are related to the dual conformal invariant integrals appearing in the calculation of the 
amplitudes. One has to look at the "parent" integral which appears while considering 
the amplitudes and shrink n propagators at the n-th order of perturbation theory. This 
dual conformal invariance together with the original conformal invariance might lead to a 
wider algebra eventually constraining the form of the answer and reveal the integrability 
property of a theory. It is important whether the powerful M = 4 covariant on-shell 
methods such as recurrent relations (see recent [21] for example) can be generalized for 
the form-factors studied in our paper. 

Note added: while finishing writing the paper we became aware of the paper which is 
closely connected to the subject studied here [47]. 
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A Feynman rules in J\f = 1 superspace 

We want to present here the essential elements of Af = 1 superspace technique relevant 
to our computations. 

In terms of M = 1 superfields the A/" = 4 SYM action can be rewritten as (hereafter 
we use the notation of [34J, see recent examples of application of the same technique in 

[351 EH E?]) 

5^= 4 = J d 8 zTr (e -flV $V v $i) + ^ J d 6 zTr(W a W a ) + ig J (fzTi ($i[$ 2 , $3]) + cc, 

(A.54) 

where the superfield strength tensor W a = D 2 (e~ gV D a e 9V ), V = V a T a is the real J\f = 1 
vector superfield and = $"T a with I = 1,2,3 are the three chiral superfields (I is 
the index of the SU(3) subgroup of SU(4)ji), T a are the generators of the gauge group 
SU(N C ) in adjoint representation. For performing SU(N C ) T-matrix manipulations we 
used FeynCalc package for Mathematica [3H]. The following normalization for T a is used 
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in which the quadratic Casimir operator 

Ti{T a T b ) = k 2 5 ab , k 2 = 1/2. (A.55) 
The relevant Feynman rules for the propagators and vertices are 

(V a V b ) = -^-6 ab ^-, 
k 2 v 

V{W<$>) = zghU^^V^j, (A.56) 

V(WV^) = 9 ^k 2 5 IJ f adm f bcm V a ^jV b ^j. 

where 5\ 2 = 5 A (9i — 9 2 ) is the Grassmannian delta function. 
The effective one-loop triple vertex is given by 

(A.57) 

where 

V = AD a D 2 D a + (p - q) a « [D a , D & ] . (A.58) 

As usual, the vertex with n chiral (anti-chiral) lines requires additional n — 1 D 2 (for 
anti-chiral D 2 ) acting on chiral (anti-chiral) lines (or n — 1 — m D 2 (for anti-chiral D 2 ) 
if m lines are external). We used SusyMath package for Mathematica [39] for performing 
.D-algebra for supergraphs. 

Traces in this case are taken over a matrices and are evaluated in D = 4 because 
dimensional reduction is used. The following set of identities is useful: 

_m ( _m\ ^.m f—m\a.B 
a = \ ) a/3 = ( ) 

Tr[l] = Tr[l] = £ (A.59) 

where D/2 = 2 in dimensional reduction and also we have 

a m a n + a n a m = -r] mn l, 

a m a n + a n a m = -r] mn l. (A.60) 
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B Scalar integrals and their e expansion 

Here we present the list of scalar integrals which we encountered in our computation 
shown in Fig. [7J All the integrals are evaluated in D = 4 — 2e dimensions. For each loop 
the factor e tlE is added in the integration measure, we also do not write 4n which always 
appear with /z 2 . 




Figure 7: The set of scalar integrals. The arc line in Gj corresponds to the presence of 
the numerator (k — p) 2 . Thick black line corresponds to off-shell leg with momentum q. 
All the other legs are on-shell. 

«o = /«^^fe^rVi + 2 + 0(e)|, (B.61) 



{2tt) d k 2 {k + p) 2 \167c 2 \fi 2 J ) \e 

i d k 1 . . 

1 = ' (2n) D k 2 (k + Pl ) 2 (k + p 2 ) 2 ( ' 62) 

16^ U J -AT 2 - 2 -^ e - 1440 e +0(e) )' 



d D k d D l 

(2n) D (2tt) d k 2 l 2 ik - Pl ) 2 {k + p 2 ) 2 {k + I - p 2 ) 2 (k + l-pi) 



^ 2 I /o_\d/o^d umii. „ \2n„ i „ \m„ i 7 „ \?n„ , ; „ \2' (B.63) 



e~ e7g Ai 2 \ V 1 / 1 5tt 2 29C 3 3tt 4 n 
16tt 2 1 /i 2 i I s\ 2 I 4e 4 24e 2 6e 32 lCj 
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(2n) D (2n) D k 2 P(k - l) 2 {k + I - Pl ) 2 (k + I + p 2 ) 2 (p 2 + I) 2 ' 

2 If 1 , 29< 3 + 4^ 





S12 V 2e 4 6e 720 
Ga = I (2n) D (2n) D k 2 l 2 {l + Pl ) 2 (k + p 2 ) 2 {k - I + p 2 ) 2 (Z - p 2 ) 2 (B-65) 

2 J_fj 7T 2 8C 3 197T 4 A 

sia 24e 2 3e 480 + (€) J ' 

d D A; rf D / 1 
5 ~ ' (2tt)^ (27r) s m\k - Z) 2 (fc - p 2 ) 2 (A; - / - Pl ) 2 (l - p 2 - Pl ) 2 ( ' 66) 

The other integrals entering into the calculations are 

/■ d D k d D l s 12 k 2 - Tr{p lP2 Ik) + Tr(pifcr P2 ) , R fi7 ^ 

5 7 (27r) D (27r)^fc 2 / 2 (fc-/) 2 (A;-p 2 ) 2 (A;-/-p 1 ) 2 (/-p 2 -p 1 ) 2 ' 1 j 

rfe ( d^J^l g 12 fc 2 - Tt(p lP2 kl) 

5 J {2ix) D {2i:) D k 2 l 2 {k-l) 2 {k-p 2 ) 2 {k-l-p l ) 2 {l-p 2 -p l ) 2 ' 1 J 

a f d D k k a @ 

Gl = J (2n) D k 2 (k + Pl ) 2 (k + p 2 ) 2 ' ^ B ' 69 ^ 

G is the scalar bubble integral, G\ is the scalar triangle integral. G 2 , G 3 G 4 can 
be computed by means of the MB representation as series in e or to all orders in e by 
means of the differential equation techniqu^f). G5 can be computed by means of the MB 
representation as series in e, the answer to all orders in e is given in [H]. G§ and G7 can be 
computed by direct evaluation of integrals over the Feynman parameters in terms of the 
hypergeometric function 2 F\ and the Appell function F\. The formulas from [42J and [45] 
were useful in verification of our computation. Gq can also be evaluated by means of the 
differential equation technique |40j, the result coincides with ours after the rearrangement 
of hypergeometric functions. Using the notation si 2 = s, su = t, S13 = u the answers for 
G 6 can be written as: 

_ r 3 (l-e)r(l + 2e) 
Cr ~ T(l - 3e) 

, d D k d D l 1 
( ' ; ' (2tt) d (2ti) d l 2 {l - k) 2 {l - Pl ) 2 (k + p 2 ) 2 (k + p 2 + p 3 ) 2 



^This integrals can be reduced to the set of master topologies presented in [40] . 
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e- 2e ~< E cr 1 1 
(16tt 2 ) 2 2?(1 -2e) 

t) 2 ' ; ^ 1 " 2e - 2 " 2e - - j> + (t) 2 ' i 1 " 2e - 2 " 2e - -t ; 

/ jj 2 \ 2e s + t + u u(s + 1 + u) 



, ^(1,1 -26,2 - 26,--^ '-)}. (B.70) 

ys + t + uy ts st 

where the hypergeometric function is given by the following expansion: (we used the Nested 
Sums computational tool jS] and HypExp package for Mathematica 



3 

2 F 1 (1 - 2e, 1 - 2e; 2 - 2e; x) = £ e n a n (x) + 0(e 4 ) (B.71) 

n=0 

log(l - ar) 
a {x) = -■ 



ai(x) 



x 

(2Li 2 (x) - (log(l - x) - 2) log(l - a;)) 



2 

a 2 (aO = flog 3 (l - x) - 31og(x) log 2 (l - x) - 31og 2 (l - x) + 7r 2 log(l - x) 

- 6(log(l - x) - l)Li 2 (x) - 6Lz 3 (l - x) - 6Li 3 {x) + 6C(3)) 

2 

a 3 (a;) = (log 4 (l— x) — 6 log(x) log 3 (l — x) — 2 log 3 (l — x) + 6 log(x) log 2 (l — x 



3x 

+ 2vr 2 log 2 (l-x) -6(log(l-x) -2)Li 2 {x) log(l-x) - 12Li 3 (x) log(l-x) 
+ 12C(3)log(l -x)- 2ix 2 log(l-a;)-12(log(l-x)-l)Li3(l-x) + 12Li 3 (x) 

+ 12Li 4 ("Trr) " 12 C(3) 
The finite part of Gq is then given by the following expression: 

,„ , 1 lit . u. « s + t + u . u(s + t + u),l 
= (lW 2 ^"^ + ^ * } " * } | • 

In the case of " 1 + 1" dimensional kinematics the following limiting expressions for G 6 
are used: 

Ge\t=o — 0, 

e" 2e7E c r 1 f 1 1 ] 

Ge\ u =o 



(16vr 2 ) 2 4e 4 { (t + u) 2 ' t 2 ' j 
e -2e 7Bcr — l l J t + s 1 t 1 1 

(16tt 2 ) 2 2e3i_2e j s (t + s) 2e ~ It 27 ~ ^ 

26 



For G 7 one has 

r f d D k d D l (k- Pl ) 2 = 

7 J {2n) D {2n) D m\i - k)\i - Pl y{k + P2 f (k+ P2 + P3 y 6 + 

- f^V) 2 ^^ 1 ' 2 ^ 1 '^^^^'^} (B.72) 

4 

2^(1, -e, 1 - e, x) = £ e" 6 n (x) + 0(e 5 ), (B.73) 



?1=0 



where 



and 



bo 


(x) = l 


h(x) ~- 


- log(l — X 


b 2 (x) 


= —Li 2 (x) 


h(x) 


= -Li 3 (x) 




= —Li 4 (x). 
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where 



2^(6,26, l + e,x) = c n (x) + 0(e 5 ), (B.74) 

n=0 
Cq(x) = 1 

ci(x) = 
c 2 {x) = +2Li 2 (x) 

2tt 2 

<%(x) = (— log(l-x)-2 log 2 (l-:r) \og(x)-A\og(l-x)Li 2 (x)-ALi 3 (l-x)-2Li 3 (x)+A( 3 ), 

c 4 (x) = -i-(47r 4 -907r 2 log(l-:r) 2 -15 log(l-a;) 4 +300 log 3 (l-x) log(x)+360 log 2 (l-a;)Li 2 (:r) 
90 

+ \og(l-x)(-3Q0C 3 +720Li 3 (l-x)+3Q0Li 3 (x))-3Q0Li 4 (l-x)-180Lu(x)-3Q0Li 4 (— X — )). 

x — 1 

In the integral G 7 there is the Appell function of the first kind defined by the following 
integral representation: 

FfahPwiV) = w ffl f duu a -\l-ur~ a -\l-ux)-^l-uy)-^ (B.75) 

1 {a)L (7 — a) Jo 
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which in our case gives us the one-parametric integral 

Fi(l,2e,l,2 + e;x,?/) = (1 + e) ^ du(l - u) e {\ - ux)- 2e {l - uy)' 1 . (B.76) 

Jo 

Expanding the integrand over e and then performing the integration one gets 

Fi(l, 2e, 1, 2 + e; x, y) = (1 + e) C du (1 + (log(l - u) - 2 log(l - ui))e 

Jo 

+ Q log 2 (l - u) - 2 log(l - ux) log(l - u) + 2 log 2 (l - e 2 ) 

+ i (log 3 (l - u) - 6 log(l - mx) log 2 (l - it) + 12 log 2 (l - «x) log(l - u) 
-81og 3 (l -ux)e 3 ) +0(e 4 ). 

Up to the second order in e it is possible to evaluate the integrals in terms of logarithms 
and polylogarithms; however, in higher orders new functions appear. 
Consider, for example, the integral 

X, = ffej^^gjM, (B.77) 
Jo 1 — uy 

where the parameters satisfy the condition < x < y < 1 . To evaluate this integral we 
use the integral representation for one of the logarithms and get 

f 1 , f 1 , — urr log 2 (l — w) _ „ . 

/ d M / da- ^ B.78 

Jo Jo (1 — uy){l — uxa) 

Now taking the integral over u one has 

2 - y-x T . -y , r 1 £*3^=i 



Ji = — log^ Li 3 -^ + 2 " ax -\ da. (B.79) 

1/ 1/ 1 - V Jo a(ax - y) 

To find the integral 

/ 7 ax ~\ da (B.80) 
Jo a(ax — y) 

it is useful to introduce a new variable 



6 = 

then 



ax 



ax — 1 



r f 85g V = -± /* , U f v do (B.81) 
Jo a(ax-y) y Jo b(l + ±=*b) K ' 



i(ax-y) yJo 6(1 + ^6) 
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and using the identity 

1 1 1-y 



6(1 + ^6) b y(l + *=*b) 

one comes to the integral 



i ( r * _ i^y , W db j (Rg2) 



l/^o 6 y ./o 1 + -^6 

The first integral is straightforward and for the second integral one can expand the poly- 
logarithm in power series and get the answer in terms of the function 

Li m , n (x,y) = Yl ( B - 83 ) 

j>i>0 J 1 

As a result one gets 

f 1 Lh^h da = _ 1 f Lk f x_\ + i / y x(l-y) \ \ ^ 

Jo a(ax — y) y\ \ 1 — xJ ' \ 1 — y 1 y(l — x) J J 

Finally putting everything together we obtain 

X, = -l iog yZl Li3 ^-l U (-^W . (B.85) 

y y 1-2/ y \ V 1— a?y \ l ~y y\ l - x )J) 

Another possibility to expand the Appell function is to use the Nested Sums computa- 
tional tool [H] which represents the Appell function as some combination of generalized 
polylogarithms. The e expansion for the Appell function then takes the form 

F 1 (l,2e,l,2 + e,x,y) = -Miz^l + I ( l og 2 (l-x) + 2 log(l-y) log(l -x) (B.86) 

y y v 

- l - log 2 (l -y)- log(l —y) — 2Li 2 {x) - 2Li 2 (|^f ) + Litivj) e 

+ - (-\ log 3 (l - y ) + -(- log(y) + 2 log(y - x) - 1) log 2 (l - y) - -tt 2 log(l - y) 
y V o 2 o 

+ log 2 (l - x)(log(a;) - \og(y) + \og(y - x) - 1) - 2Li 2 {x) - 2Li 2 f^y) 
+ log(l - x) (bg(l - y){2 - 21og(j/ - a:)) + ~ (6Li 2 (z) + 6L« 2 (\ - ~J 
+ 6Li 2 (|^f ) - 6Li 2 (y) + vr 2 )) + Li 2 (y)+2Li 3 (x)-2Li 3 ^1 - ^ -Li 3 (l-y) 
+ 2L, 3 (H) +2 ^ (S)+^s (^) -Li3(,)-C(3)) e 2 +Fin e 3 + 0(^, 
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where 

Fin = - ( -2Lii,i,i,i(l, -, 1, y) + 2 J Li 1)2jl (- ) 1, y) - 2Li 1)2 (-, y) + 2Li 3 ,i(-> I/) ( B - 87 ) 

y V 1/ 2/ y y 

+SoM+Li 3 (y)-S M-H 2!2 (y)-2Li 2X1 (l, -, y)-H 1>3 (y)+2Li lxltl (-, 1, 1, y) 

y y 

-Hi. 2 ,i{y) - S 2t2 (y) - 2Lz 1>li2 (l, -, y) + 2Lz 1)lj2 (- 1, y) + 2Lz 2)2 (- y) + # x , 2 (y) 

y y y 

+2Lz lj3 (-, y) - 5 1>3 (y) + 2Lz Mil ,(-, y) + 2Li 2tl>1 (-, 1, y) - L* 4 (y) - 2Li 2) i(-, y) 

y y y y 

+S\, 2 (y)-2Li 1 , 1 , 1 (~ l,y)+2i^^ . 

y y y J 

Here we use the following definition of the generalized Goncharov polylogarithms [30] 

Lim lr ..,m k (Xl,...,X k ) = ■^■■■^k- ( B - 88 ) 

ii>j 2 >--->«fc>0 Z l % k 

Apart from these functions we have in expansion the so-called Nielsen generalized poly- 
logarithms 

S n>p = Lin+i,i,.„,i(ar, (B.89) 

p-i 

and also the harmonic polylogarithms 

H mi ,...,m k {x) = Li mi) ^ mk (x,l,...,l). (B.90) 

fe-i 

In the case of "1 + 1" dimensional kinematics the expression for the integral Gj is 
simplified and the following limiting cases can be used: 



Gr = G7 +V + i -2e 2 J - ^ , (B.91) 



e- 2£lE cr 1 f 2 1 
(16tt 2 ) 2 2e 4 I C 

where J is the integral 

J = / rfsdy y - — r (B.92) 

Jo (tx + sy + uxy) 1+j!e 

J\t=o = 0, (B.93) 

1 r 1 tv^ 1 11 m 

J l- = -2?/„ d! '(rT^7 = -2?F7 2Fl(f - 1 + 2f ' 1 + f '-7 ) (&94) 

J U=o = -7T I / rf y7r^ — ^7 - I d y v 



2e [Jo (t + sy) 2e Jo (sy) 



2c 
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The finite part of Gj is given by the following expression: 

, „ N , _ . 1 1 f , , , u. s + u u s + u u ) 

(G 7 ) m = ( G6 , / ,„ +ITg - ¥ _ |-i + *(-) + «(-— ) - — , ^)} ■ 

(B.96) 
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